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We develop a general theory for one-way optical guiding in magnetized periodic particle arrays. 
Necessary conditions for a non-even dispersion curves are derived and presented in the context 
of Frieze symmetry-groups. It is shown, for example, that one-way guiding can be supported in 
particle strips consisting of geometrically isotropic particles arranged in transversely asymmetric 
arrays. Specific examples consist e.g. two parallel isotropic particle chains with different periods. 
The previously studied one-way effect based on the two-type rotation principle is shown to be 
a special case. In the latter the exclusion of the appropriate Frieze-symmetries is achieved in a 
single linear chain by associating a geometric rotation to each particle, thus providing the narrowest 
possible one-way waveguides. It is also shown that nearly any randomly created period may result 
in uneven dispersion and one-way guiding. 


I. INTRODUCTION 

One-way waveguides are important building 
blocks for many functional systems in acoustics 
and electromagnetics. They find applications 
as isolators and circulators, as a mean to reduce 
disorder effects in waveguides [lj, or even as a 
mean to match waveguide to a load [2] or to a 
leaky wave nano-antenna [1]. 

Non-reciprocity and one-way guiding can be 
achieved either by medium’s non-linearities, by 
spatiotemporal modulation, or by magnetic bi¬ 
asing. Non-linear non-reciprocal structures [U - 
H] use the fact that reciprocity does not im¬ 
ply symmetric field distributions of the forward 
and backward propagating modes. Therefore, 
spatially-varying material nonlinearity created 
in an initially reciprocal waveguide, may affect 
the two modes differently and eventually result 
in non-reciprocity. Since it is based on non¬ 
linear effects, this method requires large vol¬ 
umes, it is power-dependent, and may lead to 
signal distortion. Alternatively, spatiotempo¬ 
ral variation of the medium’s properties may 
violate reciprocity. For example, by using a 
traveling wave to modulate the medium’s con¬ 
stitutive parameters one may asymmetrically 
add a “momentum” (wave number) bias to 
the otherwise-symmetric mode picture. Sev¬ 
eral non-reciprocal components and one-way 
waveguides have been proposed using this idea, 
for electromagnetic waves 0,113 as well as for 
acoustic waves & These structures are linear 
and magnetic-biasing free. However, the modu¬ 
lation effect is usually quite weak thus yielding 
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a narrow one-way bandwidth and a large de¬ 
vice, many wavelengths long. The third way to 
violate reciprocity is using magneto-optical ef¬ 
fects produced by static magnetic biasing under 
which ferrites and plasmas become gyrotropic. 
Then, if a waveguide is asymmetrically loaded, 
isolation can be obtained 12]. Many of the re¬ 
cently proposed one-way schemes in optics are 
based on this idea ill Ibil - ITiit . However, these 
structures have lateral width of at least one 
wavelength. Ring resonators have been pro¬ 
posed as a mean to enhance non-reciprocity 
and reduce the device’s size. The degeneracy 
of a pair of counter propagating modes can be 
removed by azimuthal spatiotemporal modula¬ 
tion j2Ci[-[22j , or by magnetization [H] - yielding 
a mode-Q-factor dependent isolation with de¬ 
vice footprint of the order of a wavelength. 

Along an apparently different research en¬ 
deavor, linear periodic chains of equally-spaced 
identical plasmonic nano-particles were studied 
in a number of publications 01]-[Id]. Particle 
chains were proposed as guiding structures and 
unctions in [23-|32f, as surface waves couplers 
33], and as polarization-sensitive waveguides 
34j . When the inter-particle distance is suf¬ 
ficiently small, the guided optical wave in such 
chains is highly localized, i.e. it’s spatial width 
is much smaller than the wavelength that corre¬ 
sponds to the operation frequency. Hence, they 
are termed as Sub-Diffraction Chains (SDC’s). 
The chain Green’s function was calculated in 
03, and the SDC’s modal features were stud¬ 
ied using a general approach and spectral anal¬ 
ysis in [U, [§2j and were also considered in [H] . 
Green’s function theories revealing all the wave 
constituents that can be excited in these struc- 
tures, were developed and discussed in detail 
in [la]. Scattering due to structural disorder 
and its effect on the modes were studied in 03, 
and in [39]] ■ Propagation in curved SDC’s has 
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been considered recently in [4(j. The effect of 
liquid-crystal host of SDC’s has been studied 
in [4l| . where the interesting option of dynam¬ 
ically controlled SDC is suggested. 

In [ 42 U 4 H a family of magnetized plasmonic 
SDC’s one-way waveguides, based on the ‘two- 
type rotation” principle was proposed for one¬ 
way guiding and optical isolation. As has been 
shown there, exposing a conventional spherical 
particles SDC to external magnetic bias may 
formally create non-reciprocity due to the Fara¬ 
day rotation induced in the plasmonic parti¬ 
cles. However, this is not sufficient for optical 
isolation since the resulting dispersion w(/3) is 
still an even function and no preferred direction 
is created. To achieve the one-way property, 
the suggested structures consist of magnetized, 
non-spherical, plasmonic particles arranged as 
a linear array. Each particle in the array is 
rotated with respect to its neighbors, render¬ 
ing the chain chiral. The interplay between 
the field polarization rotation due to the gy- 
rotropy of the resonant particles and the struc¬ 
tural chirality is utilized in order to enhance 
non-reciprocity and create uneven dispersion. 
Significant optical isolation and one-way guid¬ 
ing are then obtained at operation points near 
the light-cone - see IT2h44I |. These structures 
have two main advantages over other one-way 
schemes. First, as said above, the effective lat¬ 
eral width of the guided mode is much smaller 
than the wavelength. Second, it has been 
shown that the required magnetic bias is con¬ 
siderably weaker than the bias used in other 
magnetization-based one-way schemes. This 
bias may be further reduced by utilizing the 
enhanced Faraday rotation in a plasmonic par¬ 
ticle consisting of a core-shell geometry such 
as the particles studied in [HI, [4(|. However, 
despite their appealing properties, a nanoscale 
fabrication of the proposed structures might be 
somewhat challenging due to the apparent re¬ 
quirement for a precisely characterized geomet¬ 
rical features in the particles. 

Here we approach the problem of enhanced 
non-reciprocity from a wider point of view. Our 
study provides the necessary conditions for un¬ 
even dispersion curves and one-way guiding in 
planar particle arrays, some examples of which 
are shown schematically in Fig. [T| The parti¬ 
cles may or may not be identical and symmet¬ 
ric. The necessary conditions for uneven disper¬ 
sion are characterized in the Frieze symmetry- 
groups framework. The previously reported 
one-way chains based on the two-type rotation 
principle [42l-l44l] are shown to be a special case 
that conveys the narrowest structure (but not 
the simplest to fabricate). One of the important 


( a ) 

A 

pi 


• • 

0000 



w x w ' 

(b) 


plml 


• • 

0000 

(c) 


pllm 


0 • 

0000 


0 0 

0000 

(d) 


piig 


• 0 

0000 


y © 

0000 

(e) 


p2 


0 0 

0 0 0 0 

O 

0 0 

0 0 0 0 

(f)| 


p2mg 


0 0 

0 0 0 0 


0 0 

0 0 0 0 

(g)J 


p2m m 


o 0 0 • • 0 0 


0 0 0 0 0 0 0 

FIG. 1. A periodic particle strip belongs to at least 
one of the seven Frieze symmetry groups, invariant 
under the following operations (in addition to in¬ 
variance under 7b): (a) No other symmetry - the 
Pi group, (b) 1Z V - the plml group, (c) TZh - the 
pllm group, (d) Glide reflection: TZhTd, d < D - 
the pllg group, (e) 7Z n - the p2 group, (f) TZ V , 
7 Z n , and TZhTd, d = D/2 - the p2mg group, (g) 
TZv, TZh , and TZ n - the p2mm group. 

conclusions of the present study is that almost 
any random construction of the structure’s pe¬ 
riod would result, under magnetic bias, in un¬ 
even dispersion and may function as an optical 
isolator. In fact, due to the structural asymme¬ 
try requirements for uneven dispersion studied 
here, and due to the fact that fabrication er¬ 
rors are unavoidable, it may turn out that un¬ 
der magnetic bias uneven dispersion is easier to 
obtain than even dispersion. 


II. FORMULATION 

We refer to the periodic particle “strips” 
shown schematically in Fig.[T| For convenience, 
let us denote the plane on which the strips lie 
as the ( x, y ) plane, x being the strip axis. The 
strip period is D. To study their symmetry 
properties we define the following translation 
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by <3, reflection about horizontal/vertical line, 
and rotation by ir operators in the (x, y ) plane 


Td{x,y) m- (x + d,y ) 

(la) 

Kh{x,y) m- (x,-y) 

(lb) 

TZ v (x,y) n- (~x,y) 

(lc) 

Kir(x,y) m- (-x,-y). 

(Id) 


We emphasize that the symmetry of any of the 
structures presented, e.g. in Fig|T] emerges not 
only from the locations of the particles, but also 
by the specific properties of the particles popu¬ 
lating these locations; e.g. their polarizabilities. 
Our aim is to study the electrodynamic prop¬ 
erties of the arrays schematized in Fig. [T| and 
their potential use for enhanced breach of time- 
reversal symmetries. In particular, we look 
for the necessary conditions under which they 
support guided modes with uneven dispersion 
curves and, eventually, one-way guiding. 

The particles are much smaller than the free- 
space wavelength A at the structure’s operation 
frequency to. Hence, we use polarizability the¬ 
ory and the discrete-dipole approximation. The 
structure’s period consists of N particles, not 
necessarily identical, each characterized by its 
own polarizability matrix a n , n = 1,..., N. In 
the present work, a necessary condition for the 
asymmetries of a n is achieved by applying a 
bias magnetic field B o = zBq , that affects es¬ 
sentially the particles dipole response in the x, y 
plane. Therefore the dipoles z components are 
ignored, rendering ct„ a set of N matrices of 
2x2 elements. We denote the location of the 
n-th particles in the m-th period by r TOi „, and 
the dipole moment excited in this particle by 
p rn n . The structure’s electrodynamics is gov¬ 
erned by the infinite difference equation 

N 

Pm,n = a nY Y^’^ r k ,e)p kti (2) 

k e=i 

where G(r,r') is the free-space Green’s func¬ 
tion matrix, by which G(r,r')p gives the elec¬ 
tric field at r due to an electric dipole p at r'. 
Here and henceforth, the absence of summation 
limits indicates a summation over all integers, 
and indicates a summation that excludes 
the v term. Due to periodicity, we have 

Pm,n = Po,ne lm/3D (3) 

hence Eq. © can be re-written as the following 
2 N x 2 N matrix equation (n = 1,..., N ) 

M(/3)P = (4a) 

N 

= (a- 1 - M£) p 0i „ - Y n M n,ePo,e = 0 
e=i 


where are N matrices of 2 x 2. The ma¬ 
trices (a- 1 — M,^) reside in the N diagonal 
blocks of 2 x 2 elements each. are 2x2 ma¬ 
trices residing in the off-diagonal blocks. 
and are given by the matrix-sums 

m£(/ 3) = ]T 0 G(ro,„,r m , n )e im ^ (4b) 

m 

M ni t(p) = Y G^o.n, r m/ ) e im ^ D . (4c) 

m 

The structure dispersion is the set of /3(w)’s 
that nullifies the determinant of Eq. (Hal) . An 
even dispersion curve is obtained if this deter¬ 
minant is invariant (up to a /3-independent mul¬ 
tiplication constant) under the map /3 H > — /3. 
Below we study the symmetries of Eq. Hal) un¬ 
der this map. 

First, we note that the matrices and 
M nA are particle-independent, and their sym¬ 
metries emerge only from the symmetries of G 
and from the symmetries (if any) of the lat¬ 
tice sites arrangements inside the strip period. 
These are studied in appendix [A] For the diag¬ 
onal blocks we have, 

m£(—/?)=m£08)= [M^(/3)] T (5) 

and for the off-diagonal blocks 

M„,,(-/3) = M,,„(/3) = [M, >n (/?)] T . (6) 

These properties hold just by a mere periodic¬ 
ity, and are proved in appendixjA] Hence, under 
the mapping /3 — /? Eq. Hal that is written 
for /3, is mapped to 

M(-/3)P= (7) 

N 

= { a n 1 - M n ) P0,n - Y " Po,t = 0 

f-1 

Hence, the mapping /3 H > — /3 is manifested by 
a block-transpose operation on the matrix that 
consists of the blocks M„^. Generally, even 
dispersion curve f3{cu) is obtained if there is, 
up to a multiplication constant, a determinant¬ 
preserving transformation by which one can ob¬ 
tain Eq. © from Eq. Hal) . Below we explore 
the conditions for even dispersion. The simulta¬ 
neous violation of all these conditions is a nec¬ 
essary condition for uneven dispersion. 

A. Symmetric a n ’s 

If the a n ’s are symmetric 

at n = a^, V n (8) 

then, from Eq. ([6]) we have [M(/3)] T = M(—/?). 
Hence |M(/?)| = |M(— (3)\ (transpose preserves 
the determinant), leading to even dispersion. 
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B. Non-symmetric and different a„’s 

Assume now that a n are non-symmetric and 
also not necessarily identical. For simplicity, let 
us number the polarizabilities within a period 
along the positive x and y directions; i.e. from 
left to right and from bottom to top. With 
the use of the inversion matrix I 1 2 at defined in 
Appendix lBl Eg. (IB2I) . we have 

M(/3) = (9a) 

if the following conditions are satisfied: 

a-N+i-n = a n , M® +1 _ n = , Vn (9b) 

and 

M£ >ra = MAr + i_ raj jv + i_£, V£, n. (9c) 

These conditions imply inversion symmetry (in¬ 
variance under TZ n ), that applies on both the 
particle coordinates [see appendijAl Eqs. (IA8I) 
m, as well as on the particle properties. 
Since I* 2 jv is determinant preserving, it follows 
that in inversion symmetric chain structures 
the dispersion curves are always even, indepen¬ 
dently of the asymmetric form of each of the 
individual polarizabilities. 

To find another condition for even dispersion, 
we refer to Eq. (IA4bD . by which applying 7 Zh on 
the chain is equivalent to matrix transformation 
by I h 2 . Consequently, from Eqs. (l4bl) (l4cl) it 
follows that applying IZh on the structure is 
equivalent to the matrix transformations 

M^I h 2 Mfl h 2 (10a) 

M/, n £* I h 2 M e ,n I h 2 (10b) 

this operation merely changes the signs of the 

two off-diagonal terms in the matrices and 

M £ tn . We refer now to Eqs. (lAllal) (lAllbl) 
that hold if the set of points in a period pos¬ 
sesses horizontal reflection symmetry. Hence, 
using the horizontal inversion matrix operator 
I h 2 JV defined in appendix [Bl we have 

M r (/3) = I h 2ArM(-/3)I h 2iv (11a) 

if the structure possesses horizontal inversion 
symmetry and if 

a^ = I h 2 a„I h 2 , Vn = 1,... TV. (lib) 

Since the operations in Eq. (Illal) are deter¬ 
minant preserving, it follows that under these 
conditions the dispersion curve is even. Note 
that the condition in Eq. (lllbl) applies not only 
to conventional spherical particles, but also to 
plasmonic particles under magnetization. 


Finally, we note that the glide-reflection 
transformation is defined as TZhTd - see Fig. [TJl. 
However, Td has no effect on G and subse¬ 
quently no effect on and M^„. There¬ 
fore the last results apply to symmetries under 
glide-reflection transformations as well. 

To summarize, it has been shown that if 
the structure possesses inversion symmetry (p2 
group), reflection symmetry with respect to the 
chain axis (pllm group), or glide-reflection 
symmetry (pllg group), an even dispersion is 
guaranteed even if the particles themselves pos¬ 
sess non reciprocal property such as that arises 
e.g. by magnetized plasmonic material. 


C. The non-uniqueness of the unit cell 

So far we have discussed the symmetries on 
the unit cell level. However, while the period 
length, D, is unique, the choice of a unit cell 
is not. The latter always consists of N parti¬ 
cles, and there are up to N different choices for 
the unit cell and even more possible M (/?) (due 
to the freedom of ” numbering” the particles in 
the same unit cell 0)- Let us define the set 
M as the set of all possible matrices M(/3) for 
a certain structure. If the determinant of the 
matrix M((3) vanishes for a certain choice of 
the unit cell, then it must vanish for all M - 
since all the matrices in M represent the wave 
dynamics of the same structure. Therefore, if 
the operations done on M ((3) in the previous 
subsections do not result in the mentioned vari¬ 
ations of the same M ((3) but do result in these 
variations of one of the other matrices in M 
then the structure must still possess even dis¬ 
persion. For example, if the unit cell does not 
posses an inversion symmetry, but upon invert¬ 
ing the unit cell we obtain a differently struc¬ 
tured unit cell, but of the exact same original, 
pre-inverted structure - then the dispersion is 
rendered even since the substitution (3 —> —{3 
still preserves the zero determinant. A good 
example for this is the pllg group, displayed 
in Fig. [1]1. A unit cell does not posses any of the 
discussed symmetries - invariance under TZ^ or 
under 7 Zh- However, under 7 Zh we obtain a dif¬ 
ferent unit cell of the same structure, and there¬ 
fore it is reciprocal. This fact helps us finally 
establish the list of Frieze groups that host even 
dispersions: pllm, pllg, p2, p2mg, p2mm. 
Uneven dispersion can only be found in groups 
Pi and plml. 
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D. Random period 

Consider a periodic chain whose period is 
randomly generated (e.g. randomly locating N 
particles inside the period). Clearly, for N = 1 
the notion of randomness becomes meaningless. 
For N = 2 and if aq = a 2 , one may verify 
that the period would always be invariant un¬ 
der lZ n , i.e. it resides in the p2 group and hence 
it possesses an even dispersion. Note that this 
argument holds only for dipole-moment excita¬ 
tion. If higher order excitations become rele¬ 
vant (e.g. quadrupoles), then the above conclu¬ 
sion may not hold. However, for N > 3 nearly 
any random construction would reside in the 
Pi group whether the particles are identical or 
not, yielding uneven dispersion under magnetic 
bias (see Sec. Mil belowh 
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FIG. 2. Examples of symmetry considerations 
in isotropic and non-isotropic particles. (a) A 
structure made of non-isotropic particles that have 
isotropic (scalar) polarizability when not mag¬ 
netized. (b) Non-isotropic particles with non¬ 
isotropic polarizabilty. 


E. Nonreciprocal sector-way metasurfaces 

A new family of metasurfaces, termed as 
metaweaves, can be obtained by invoking a pro¬ 
cedure of “weaving” one-way particle chains 
(one-way “threads”). The resulting surface 
possesses a novel “(j> sector-way” dynamics: 
when the surface is excited by a point dipole it 
allows the propagation of trapped modes only 
into a cone whose vertex-angle is </>. This struc¬ 
ture has been suggested and studied in [50l ]. 
With the present approach, however, one may 
design non-reciprocal sector-way metasurfaces 
without the need to directly “weave” one-way 
threads. Generally, the symmetries of such 
metasurfaces can be studied in the framework 
of the wallpaper groups in a manner similar to 
the study in the previous subsections. However, 
there are 17 symmetry groups, hence for brevity 
we use here a rule of thumb (see example be¬ 
low). To achieve sector-way dynamics it is suf¬ 
ficient to create a periodic metasurface whose 
period does not possess IZh, 7Z V , and 1Z„ sym¬ 
metries (i.e. the two-dimensional period should 
not be invariant under these operations). 


rotated by an angle smaller than 7r/4). How¬ 
ever, since the polarizabilities are isotropic, this 
asymmetry will not show up in the particle- 
dipole model, and the corresponding dispersion 
is even also under magnetization. It should be 
emphasized that the theoretical analysis above 
relate only to the particle’s polarizabilities, as 
seen in e.g. Eq. (J5J) , Eq. (l9bl) . and Eq. (Illbl) . 
Hence, the symmetries demonstrated in Fig. [T] 
should be perceived as symmetries that relate 
to particle’s locations and polarizabilities (and 
not particle geometries). In light of the above, 
the structure in Fig. [2ji would possess uneven 
dispersion only if the particles are sufficiently 
large, so that higher order multipoles are sig¬ 
nificantly excited and are needed to model its 
dynamics. To contrast, the structure in Fig. [2}o 
consists of particles that possess anisotropic po¬ 
larizability. Hence, it resides in pi group both 
in the formal geometrical sense (like the cubes 
structure), and in the polarizability sense (un¬ 
like the cubes structure). 


F. Distinctive features related to 
modeling of non-reciprocal particle arrays 

Consider a particle with no circular symme¬ 
try, but with an isotropic (scalar) polarizability 
when not magnetized. A cube made of isotropic 
dielectric material is an example for such a par¬ 
ticle HH. An example of a structure made of 
such particles is shown in Fig. [2Ji. Formally, 
this structure resides in the plml group (or in 
Pi group if the particles in the a lower chain are 


In conclusion of the above discussion, par¬ 
ticle structures may possess “hidden” symme¬ 
tries that take effect only when the particles are 
small, and disappear when considering larger 
particles or continuous medium. Furthermore, 
the ability to assign “point-wise” rotation to 
a particle whose polarizability is non-isotropic 
(e.g. ellipsoids), adds degrees of freedom not 
available in continuous medium or in ID pho¬ 
tonic structure such as the one studied e.g. in 
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III. EXAMPLES 


A. plml group 


All the structures in the examples below con¬ 
sist of plasmonic particles, exposed to a uni¬ 
form transverse magnetization bias B 0 = zB 0 . 
Generally, the single particle polarizability un¬ 
der magnetization can be found, e.g. in [43[ 
for metal particles of ellipsoidal shape, and 
in [H] for engineered Graphene particles us¬ 
ing more detailed quantum model based on 
the Kubo theory [49|. We note that the lat¬ 
ter possesses quasi-static resonance very simi¬ 
lar to that of the Drude model, with resonant 
frequency whose role is similar to that of ui p in 
metals, but in the 40-50 THz frequencies. For 
all particles the off-diagonal terms of a n are 
imaginary and of opposite signs, i.e. they satisfy 
Eq. (lllbl) . The strength of the non-reciprocity 
in a is measured generally by the parame¬ 
ter k = \a. x y/a xx \. Furthermore, note that 
k = tand where 9 is the polarization-rotation 
angle associated with the non-reciprocity. For 
metal particles under the Drude model n gets 
a simple algebraic expression, from which the 
strength of the magnetization level is easily ex¬ 
tracted: k = ujb/ojp where co p is the plasma 
resonance frequency and u>b = eBo/m e is the 
cyclotron frequency. For the engineered par¬ 
ticles in [45j the dependence of n on B 0 is 
much more complicated - see e.g. Eqs. (10)-(11) 
in the supplementary information of |45j . To 
get a feeling of the numbers involved, we note 
that in the examples below we typically have 
k ~ 0.0075 — 0.01. Under the Drude model and 
for metals such as Ag, this yields B 0 way be¬ 
yond any practical realization. However, in the 
particles of [45[ , the magnetization strengths re¬ 
quired for these values of k is significantly less 
than 1 Tesla in the relevant frequency regime 
(see e.g. Figs. 2,3 there). 

To keep the algebra simple, the polarizabil¬ 
ities of the magnetized spherical and elliptical 
particles in all the examples below were calcu¬ 
lated according to the formulas provided in [431 . 
However, the final algebraic form of a in [451 ] is 
precisely the same and it needs Bq < 1 Tesla. 

Finally, we note that the structures in the 
present study are inherently non-Bravais lat¬ 
tices. Hence, they possess 2N dispersion 
branches for each plasma resonance; N is the 
number of particles in a period, and the factor 
2 is due to the vector-nature of the particle’s 
dipole moment. To avoid cluttering the disper¬ 
sion diagrams we show below only the branches 
where the non-even dispersion was examined 
and one-way guiding was obtained. 


1. By virtue of particle geometry 

Figure [3[i shows a structure consisting of 
two parallel chains of plasmonic particles, with 
identical period lengths. The chains differ by 
the corresponding particle’s volume: the par¬ 
ticles in one chain (lower) are of volume V, 
while those of the other chain (upper) are of 
volume £V, £ > 0. The corresponding polar¬ 
izabilities scale accordingly (a is proportional 
to V). Clearly, for £ ^ 1 the structure be¬ 
longs to the plml group. The other param¬ 
eters are D = A p /3, a = A p /12, do = 0.8 D, 
where X p is the plasma resonance wavelength 
and a is the particle diameter. The magnetic 
bias strength corresponds to k = 0.01 ( =bOb/oJ p 
in Drude metals). The dispersion of plasmonic 
modes guided by the structure is obtained by 
looking for the real values of that nullify 
the determinant of Eq. dial) . Figure [3):> shows 
the dispersion for £ = 0.2 (the particle’s di¬ 
ameter scales as 0.2 1 / 3 « 0.585). The uneven 
nature is evident. To get a feeling of how it 
depends on the parameter £, Fig. [3J: compares 
/3(wi) on the left side (/ 3l ) and on the right 
side (/3r) of the dispersion, at a frequency 
for which both are real. Only at £ = 1 one 
obtains (3l = Pr- Finally, Fig. 0:1 shows the re¬ 
sponse of the structure when excited by a point 
dipole located at its center, with the frequency 
shown by the horizontal line in Fig. [3J:>. Note 
that at this frequency there is only one inter¬ 
section with a guided mode. This mode has 
negative /? and negative group velocity (the red 
line in Fig. [3 Jd). There are also two intersec¬ 
tions with the light-line cone corresponding to 
the so called “light-line modes”. The latter are 
practically unexcitable in SDC-type structures, 
whether reciprocal [Hi [38[, or non-reciprocal 
[HI . Hence only —x propagating mode can be 
excited, as seen in the response in Fig. [3ji. 


2. By virtue of lattice geometry 

Figure IIJj, shows another structure of the 
plml group, consisting of two parallel chains 
of identical particles but with period lengths D 
and D/ 2. The other parameters are D = A p /4, 
a = A p /24, do = 0.5 D, and k = 0.0075 (=Wf,/w p 
in Drude metals). The uneven dispersion and 
one-way operation are seen in Figs. |3)d,c. 
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FIG. 3. A planar periodic structure of group plml 
possessing uneven dispersion under magnetic bias 
normal to the structure’s plane, (a) The structure 
geometry, (b) The dispersion for £ = 0.2. The 
vertical blue lines are the light cone, (c) A mea¬ 
sure of the uneven dispersion, (d) The response 
to a point dipole located at the center, with the 
frequency shown by the horizontal blue line in (b). 
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FIG. 4. (a) Two parallel chains of identical parti¬ 
cles, but with different periods, (b) The structure 
dispersion (c) The response to a point dipole lo¬ 
cated at the structure’s center, with the frequency 
shown by the horizontal blue line in (b). 


B. pi group 

In this group the structure period possesses 
no symmetry at all. Hence, in principle, it is 
the easiest to fabricat or synthesize. In fact, 
nearly any period that consists of randomly lo¬ 
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FIG. 5. A planar periodic structure of group plml 
if 5 = ±0.573, and to the pi group if 8 ^ ±0.573, 0. 
For 5 ^ 0 it possesses uneven dispersion under mag¬ 
netic bias normal to the structure’s plane, (a) The 
structure geometry, (b) The structure’s dispersion 
for S = 0.4d. (c) A measure of the uneven disper¬ 
sion. (d) The response to a point dipole located at 
the structure’s center, with the frequency shown by 
the horizontal blue line in (b). 


cated particles , would result in uneven disper¬ 
sion under transverse magnetic bias, and conse¬ 
quently in one-way operation. Below we show 
two examples of ordered and random periods. 


1. By virtue of lattice geometry 

The geometry shown in Fig. [5jr consists of 
three parallel chains of identical particles and 
equal period length 73. Two of which are pre¬ 
cisely aligned, while the third is shifted along 
its axis by the distance 6. Note that for 8/D = 
±0.5 the structure belongs to the plml group, 
and for 5/D ^ ±0.5,0 the structure is in the 
Pi group. The parameters are: period length 
D = Ap/6, particle radius a = A p /24, distance 
between the chains do = 0.7573, and magnetic 
bias level that corresponds to k = 0.01. Figure 
[SJ) shows the dispersion for d = 0.473. Figure 
Et shows \Pl\ — |/3r| at a frequency for which 
guiding exists for both sides (wi = 0.556053w p ), 
vs. S/D for the range —0.5 < S/D < 0.5. Fig¬ 
ure [SJl shows the structure’s response to a point 
dipole located at its center. 

To get a feeling of actual array parameters 
and the effect of loss, we note that when us¬ 
ing Ag particles the parameters above corre- 






































FIG. 6. The structure in Fig. [5] but for Ag particles 
with loss, (a) Response, (b) The structure isolation 
ratio \p{— *)| / |p(*)| that provides a measure of the 
one-way property. 

spond to A p = 150nm [HTJ] , particle diameter 
is 2 a = 12nm, D = 25nm, and do = 19nm. 
The loss parameter of Ag is r = 0.5 x 10 -12 
Sec. Due to the presence of loss /? becomes 
complex and dispersion curves such as those in 
Fig. [5] cannot be drawn. Instead, we use the 
working point shown in Fig. [5] to simulate the 
structure response. The results are shown in 
Fig- El It is seen that while the signal itself de¬ 
cays as it propagates along the structure, the 
one way property is still profound and survives 
the loss, as evident from the the isolation ra¬ 
tio. Achieving the required k under the Drude 
model for Ag requires magnetization levels that 
are far beyond practical values. However, we 
note that essentially the same physical dimen¬ 
sions would apply for structures made of the 
particles studied in [45[ (2a ~ 10 to 20 nm and 
central frequency in the range of 40-50 THz), 
with Bq in the order of 1 Tesla or less. 


2. By virtue of particle geometry and randomness 

In all the examples so far, our structures 
consist of at least two parallel periodic chains. 
Can one obtain uneven dispersion using a single 
chain? It is clear that as long as the particles 
themselves are spherical, invariance of a single 
perfectly linear chain under 7 Zh cannot be bro¬ 
ken. Hence such structures belong to the pllm 
group (at least) and their dispersion would al¬ 
ways be even. 

However, if the particles are not precise 
spheres (e.g. ellipsoids), one may assign the 
property of geometrical rotation to each parti¬ 
cle. This property is not invariant under reflec¬ 


tion, hence it opens the way to achieve uneven 
dispersion using a single chain. This approach 
of using chiral chains to achieve one-way plas- 
monic guiding has been suggested and studied 
in |12H l 41. Since it incorporates geometrical ro¬ 
tation and Faraday rotation simultaneously, it 
is termed as the tw o-typ e rotation principle. We 
refer the reader to [4J, [44j for examples of pla¬ 
nar chains. Hence, in the context of the present 
study, the two-type rotation principle can be 
viewed as a method to generate pi or plml 
structures systematically, using only a single 
chain. The studies in [43, [44| dealt with struc¬ 
tures where the rotation of the n-th ellipsoidal 
particle is given by 9 n = nA9. If A9/t: is ra¬ 
tional, this would result in a perfectly periodic 
chain. 

The generalization of the previous studies, 
as suggested by the present work, is not only 
of formal nature. It may have also impor¬ 
tant practical implications. In fact, almost 
every random generation of a period, would 
end up in pi group and may provide one-way 
guiding. We demonstrate this important find¬ 
ing within the family of the two-type rotation 
structures. Figure [7}i shows, schematically, a 
chain of rotated prolate ellipsoids whose pe¬ 
riod D consists of three particles. The ellip¬ 
soids rotation angles are chosen at random. In 
the specific example here we have ($i, # 2 , 9 3 ) = 
(0.329, 0.794, 0.468)7r. Other parameters are: 
D = Ap/3, a = \p/27 (prolate’s major axis), 
b = 0.9a (prolate’s minor axis), and k = 0.0075. 
Figure [TJ? shows the corresponding dispersion, 
and Fig. [TJ 2 the chain response to a point dipole 
excitation at its center. One way behavior is 
evident, with isolation of nearly 10 -4 over dis¬ 
tances of 5-6 wavelengths. 


C. Nonreciprocal sector-way metasurfaces 

An example of a sector-way metasurface, or 
meta-weave, is shown in Fig. [HJi. It is a gen¬ 
eralization of the structures suggested in [50| . 
The structure period in both directions is D, 
and it consists of only three spherical parti¬ 
cles located at rq = (0,0), r 2 = (0.75,0)11, 
and r 3 = (0.25,0.75)11. Clearly, this period 
does not possess TZh,TZ v , and TZ^ symmetries. 
We have calculated the dispersion surfaces and 
dipole response of the structure with the fol¬ 
lowing parameters: D = A p /2, a = A p /24, and 
k = 0.005 (the weakest magnetic bias in our 
examples). Figure |8 }d shows the dispersion con¬ 
tour for the frequency oj = 0.5762249w p . Note 
the asymmetric shape: it is not invariant under 
the map (3 1 —t —/3. The circle in the center is the 
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FIG. 7. An example of a periodic chain of ran¬ 
domly rotated ellipsoidal particles, (a) A schematic 
view of the structure. (b) The dispersion for 
{ 61 , 62 , 63 ) = (0.329, 0.794, 0.468)tt. (c) Response 
to a dipole excitation. 
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FIG. 8 . The metaweave, (a) Lattice geometry, pos¬ 
sessing no symmetry under 1Zh,R-v, and TZ n . Un¬ 
der normal magnetic bias this structure possesses 
sector-way guiding, (b) A dispersion contour. The 
red circle in the center is the light-cone, (c) Re¬ 
sponse to a point-dipole excitation. 


light-cone. Figure [ 8 J 1 shows the response to an 
excitation by point dipole located at the center, 
with the frequency corresponding to Fig. [5 }d. 
Sector-way response is evident. 


General symmetry considerations were used 
to predict the conditions under which one may 
obtain uneven dispersion curves an optical iso¬ 
lation in particle arrays. The idiosyncracies of 
dipole-dipole interactions, related to asymme¬ 
try in particle arrays were discussed. It has 
been shown that under magnetic bias, struc¬ 
tures belonging to the pi and plml Frieze 
groups shown in Fig. [1] possess this disper¬ 
sion asymmetry. The uneven dispersion can 
then be utilized for optical isolation. The sin¬ 
gle chain structure for one-way guiding, based 
on the two-type rotation principle, is shown to 
be a special case where one excludes the neces¬ 
sary symmetries by assigning geometrical rota¬ 
tion to each lattice point. Hence, it provides 
the narrowest optical isolator. It is further 
shown that almost any random periods with 
more than two particles would result in uneven 
dispersion. 

The results of this work provide a set of rules 
that can be useful for designing optica isola¬ 
tors, and as building blocks for optical circu¬ 
lators. Since one-way guiding inherently su- 
presses propagation in the “wrong” direction, 
applications that require minimization of back- 
reflections may benefit. These applications can 
be, for example, new ways to feed electromag¬ 
netic components as antennas and resonators in 
a perfectly matched fashion 09- 


Appendix A: Matrices symmetry properties 


The free-space electric field E{r) due to a 
point dipole p at r' in is given by 

E(r) = Eq 1 G{R) |/c 2 (n x p) x n 

+ [3n(n • p) - p] (Jp - | (Al) 

where G{R) = ( 47 t R)~ 1 e lkR , R = \r — r'\, and 
n = ( 1 — r')/R. From the above, the free-space 
Green’s function matrix used in e.g. Eq. m can 
be written as 


G(r,r') = e^G{R) 


k 2 (I - U) 




(A2) 


where I is the identity matrix, and U = h T n. 
G(r, r') satisfies the following symmetry rela¬ 
tions 


G(r’, r') = G T (r, r') = G(r', r) (A3) 
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that are a manifestation of the free-space reci¬ 
procity and symmetry. If r, r' are restricted to 
the z = 0 plane, then G is reduced to a sym¬ 
metric 2x2 matrix 

G{r,r')=( 9xx <J p> ) , (A4a) 

\ yxy yyy J 

satisfying 

G (lZ h r,TZ h r') =( ~ 9xy ) 

\ yxy yyy J 

= I h 2 G(r,r')I h 2 (A4b) 

where I h 2 = diag(—1,1). The last equality 
above deals with reflection about the x axis. 
Similar relation holds for the operation 1Z V . 
Furthermore, inversion leaves G unchanged 

G(7^r,7^r') = G(r,r-') (A4c) 

We turn to discuss the matrices and 
M„^, defined in Eqs. (libl) (l4cl) . Due to the 
periodicity of our structure we have, Vm,n,£ : 

r 0 ,n - r m>n = r_ TO)< - r 0 ,e = -xmD (A5) 

The first symmetries of and M n> £ emerge 
only from Eq. (IA3I) and from the strip period¬ 
icity. For M n(0) we have, 

M,?(-/3) = y°G(r 0 , B ,r m , Tl )e-^ fl 

m 

= ^°G(ro,n,r_ m ,„)e im/3i? (A6) 


r*o,i, f*o, 2 j • • •, T*o,iv possesses an inversion sym¬ 
metry if Vn = 1,..., iV, 

0 ,n = U0,n) — ^0,iV-|-l — n- (A8) 

Using this result together with Eq. (IASI) . 
Eq. (IA4cl) . and Eqs. (l4bl) Bel) , we obtain 
for inversion symmetric set of points ( n , £ = 
1 

M^ +1 _„ = (A9a) 

M i tn = M7v + i_ raj jv + i_^. (A9b) 

Likewise, if the set of points possesses horizon¬ 
tal reflection symmetry 

Rh^0,n = (^0,n; 2/0,n) = 'l'0,N+l—n (AlO) 

then, using Eq. (IA4bl) in the definitions of 
and M„ £ we obtain 

I h 2 M^I h 2 = M£ +1 _ n (Alla) 

I h 2 M„yI h 2 = Mjv+i-n.jv+l-/ (Allb) 

Appendix B: Inversion and reflection 
operations 

Let R presents a set of N points 
in the plane; R = (ri, r 2 ,... rjv) = 
(■ xi,yi,X 2 ,y 2 ,---,x n ,y N )• Assume that 

the points are numbered, e.g. from left to 
right and from bottom to top. Then, up 
to an arbitrary linear shift of the entire set, 
an inversion TZ^ of the set of points can be 
characterized by the matrix I* 2 jv 

TZ-ttR = I 1 2nR T (Bl) 


We use now Eq. (1A5D with £ = n, and Eq. (1A3I) 
in Eq. (lAfil) . The result is Eq. (|5]i. 

Regarding M n> <?(-/3), 


however, from Eq. 

ro t £. By using this result and Eq. 

Eq. JA7|, we obtain Eq. ©■ 

Inversion and reflection symmetries If 
the structure possesses more symmetries than 
that of a mere periodicity, their footprint ap¬ 
pear as more symmetries of M^ 1 and M„ j. To 
explore them let us number the points within 
a period from left to right and from bottom 
to top. Then, the set of points in a period 


where the non-zero entries of P 2 jv occupy only 
the secondary block-diagonal elements, and are 
given by N conventional 2x2 identity matrices 
I 2 [I 2 = diag(l, 1)], 


(B2) 


We note that |P 2 jv| = 1, hence operation by 
I 1 is determinant preserving. Likewise, we note 
that reflection about the x axis-horizontal re¬ 
flection 7 Zh is characterized by the matrix I h 2 jv 




(0 0 . 

. 0 

h\ 

A7) 


0 0 . 

I 2 

0 

I"2iV = — 




0 

0 



\h 0 . 

. 0 

0 / 

T—m.l — Tm,n 






lZ h R = l h 2NR J 


(B3) 


where I h 2 7 V = diag(I h 2 , I h 2 ,..., I h 2 ), and 
I h 2 = diag(—1,1). 
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